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ABSTRACT 

The primary purpose of this paper is to provide a real- 
istic estimate of the errors introduced in microwave rang- 
ing measurements as a result of diurnal, seasonal, and 
annual changes in the earth T s atmosphere. This analysis, 
which is based upon the extensive atmospheric investiga- 
tions described in the literature, is directly applicable to 
such systems as the Unified S-band System which will be 
employed in the forthcoming Apollo Missions. The effects 
of both the troposphere and ionosphere upon measurements 
of range, range- rate and elevation angle are discussed. 
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ATMOSPHERIC TRACKING ERRORS AT 
S> AND C-BAND FREQUENCIES 


by 

P. E. Schmid 

Goddard Space Flight Center 


INTRODUCTION 

During the past decade considerable effort by many independent investigators has been directed 
toward the theoretical study and analysis of the effect of the atmosphere upon microwave radio 
propagation. A number of mathematical models of the atmosphere have been devised in an attempt 
to estimate the influence of the earth T s atmosphere upon measurements of range, range rate, and 
elevation angle of earth satellite and space vehicle tracking systems. To a certain extent the sys- 
matic or biasing effects of the atmosphere are predictable and can be corrected for if sufficient 
data regarding the atmospheric state are available. This paper examines the error, attributed to 
an imperfect knowledge of the atmosphere, that occurs while using microwave energy to observe 
spacecraft at altitudes greater than 200 km above the earth. To the extent possible experimental 
verifications of the atmospheric models cited are presented. 

The effects of the troposphere and ionosphere upon radiowave propagation can best be predicted 
if an accurate refractivity profile of the atmosphere is available. Thus one approach to obtaining 
correction data is based upon atmospheric refraction profile measurements just prior to a tracking 
operation and a subsequent ray tracing by means of a digital computer (Reference 1). 

A second approach is to base tropospheric corrections on surface refractivity measurements 
(which can again be performed just prior to tracking) and, as discussed in the section entitled 
’’Ionospheric Errors,” upon the average state of the ionosphere. In this approach a mathematical 
model of the atmosphere must be used and the correction is, of course, no better than the model. 
The accuracy of the National Bureau of Standards Central Radio Propagation Laboratory exponential 
model troposphere has been verified by means of experiments which permit a comparison of cal- 
culated and measured angles of arrival of solar microwave radio energy. 

Figure 1 shows the mean total measured refraction of solar radio energy at 16.2 kMc and 
34.5 kMc compared with that calculated with the exponential model. The 16.2 kMc measurements 
(References 2 and 3) were conducted to examine the accuracy of the National Bureau of Standards 
exponential model prediction equations through a comparison with some precise measurements of 
total absolute refraction made with a radio sextant by the Collins Radio Company. This radio 
sextant consisted of a high-gain, narrow -beamwidth antenna which automatically tracked the radio 
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noise emanating from the sun's disc by means of 
continuous circular scanning. This principle is 
a common technique in radar automatic angle 
tracking whereby the antenna feedpoint is mu- 
tated (either mechanically or electronically) to 
cause a slight excursion of the spacial antenna 
pointing while the main antenna structure or dish, 
in the absence of any received signal at the an- 
tenna aperture, remains fixed. In the 2 cm wave- 
length microwave region, the sun radiates ther- 
mal noise energy similar to a blackbody having 
a temperature of approximately 7000 °K. With 
the output of the high-gain antenna connected 
to a sensitive microwave receiver, the maxi- 
mum signal output occurred when the scanning 
antenna pattern was centered upon the sun's 
disc. As the antenna scanned off axis, the re- 
ceiver output dropped to zero and an error 
voltage was generated which caused the antenna 
drive motors to redirect the antenna beam 
toward the sun’s center. 

The assumption that radio energy is uniformly generated over the sun's surface is valid except 
during times of anomalous sunspot activity. Fortunately such sunspot activity is easily discerned 
by the nature of the receiver output and consequently any data recorded during such periods was 
disregarded. The radio sextant- measured pointing was compared to calculated positions of the sun 
based on highly accurate ephemeris data supplied by the U.S. Naval Observatory. A thorough 
discussion of the error analysis associated with the approximately 40000 separate refraction meas- 
made over a five month period with the radio sextant is presented in Reference 2. A thorough dis- 
cussion of the exponential mathematical model for the troposphere as well as related topics is 
given in Reference 4. 

A third approach to atmospheric correction is to use average atmospheric data without regard 
for diurnal or seasonal meteorological fluctuations. The latter technique, as will be shown, can 
result in appreciable measurement error at low elevation angles. 

The angle of arrival measurement of radio energy from the sun is an example of a test de- 
signed to investigate the absolute accuracy of the exponential mathematical model of the troposphere. 
The refraction in this case is attributed solely to the troposphere since the solar radio energy at 
16.2 kMc is well above those frequencies affected by the ionosphere. The exponential model can 
also be checked by comparing calculated and measured refractivity profiles. Figure 2 shows a 
typical measured profile obtained at Valkaria, Florida, by means of a radiosonde (Reference 5) 
compared to the exponential model based on a measure of the refractivity at the surface and at a 
height of 1 km above the surface. The two curves are seen to coincide from 0 to 5 km and show a 
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0 MEASURED REFRACTION OF 16.2 kMc 

SOLAR RADIO ENERGY, 1959 
Nj =332 (Reference 3) 

V — MEASURED REFRACTION OF 34.5kMc 

SOLAR RADIO ENERGY, 1956 (Reference 13) 

3— CALCULATED REFRACTION 

EXPONENTIAL MODEL OF ATMOSPHERE 
h > 450km N s =313 
RAY TRACING (Reference 4) 
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Figure 1 —Atmospheric refraction of solar radio energy. 
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maximum deviation, AN, of about 10 at 12 km. 

Field tests have indicated that refraction cor- 
rections based upon refractometer profiles, 
radiosonde profiles and exponential model pro- 
files are calculable to about 1.6 percent on a 
statistical basis (Reference 6). In terms of ele- 
vation angle error remaining after attempted 
correction, this amounts to ±0.04 mr at 5° ele- 
vation. In the section entitled "Tropospheric 
Errors" it is shown that this error decreases 
rapidly with increasing elevation. The statisti- 
cal fluctuations in the correction procedure are 
largely due to the time-dependent variations of 
the refractive index profile. The variations in 
surface refractivity are directly related to 
overall tropospheric profile variations, and 
hence to the error introduced in ranging meas- 
urements, by means of the exponential profile 
ray tracings of B. R. Bean and G. D. Thayer (Reference 4). 

In most of the tropospheric models described in the literature, the effects upon radio wave 
propagation are related to surface refractivity, N s . Thus, if the model is appropriate at a particular 
value of N s and a compensating correction is made to ranging data, it is of interest to note the de- 
viations from this correction, which might be anticipated as N s , slowly fluctuates throughout the 
year. The variations come about as the changing seasons alter the air temperature, total air pres- 
sure and air water vapor content. The same parameters are also affected over a much shorter 
cycle as the air mass is subjected to daily temperature and moisture variations, especially those 
associated with the transition from local day to night. Thus as shown in Figure 3, the overall 
pattern of refractivity change during a span of a year is one of slow seasonal variation modified by 
more rapid diurnal cycles. In the United States the months of May through August are generally 
those of greatest refractivity fluctuation and the predominant change in N s during any given month 
at a particular site can be attributed to diurnal effects. 

The National Bureau of Standards has com- 
piled a book of surface refractivity charts and 
data based upon 8 years of surface weather data 
from 60 weather stations in the United States 
and 5 years of surface weather data from 306 
worldwide stations (Reference 7). The surface 
refractivity values tabulated in this publication 
are derived from measured air temperature, 
total air pressure and water vapor content. The 
basic tabulated data for each site within the 
United States is the 8 year refractivity average 



Figure 3—Typical variation in surface refractivity (data 
from Reference 7, for Miami, Florida). 
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Figure 2— Typical measured tropospheric index of 
refraction profile compared to exponential model. 
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for a particular time of day during a given month. For example, the 8 year mean value N s at 
Miami, Florida, for 12 noon during the month of August was 369.4 with a standard deviation of 8.9. 
This distribution is based on 248 separate refractivity calculations. The mean value as calculated 
for the even hours of the day clearly traces out the average diurnal cycle associated with each day 
of a particular month. This average diurnal variation is indicated in Table 1 for seven weather 
stations which include areas near the tracking stations at Cape Kennedy, Florida; White Sands, New 
Mexico; Corpus Christi, Texas; Blossom Point, Maryland; Goldstone, California; and E. Grand 
Forks, Minnesota. 

The yearly variation in N s for each site is presented as the difference between the maximum 
and minimum monthly value of N s . For example, for Miami, Florida, the maximum monthly mean 
is 381.0 (September for 2200 hours) and the minimum is 336.5 (January for 1400 hours); there is 
a yearly average change of 44.5. These variations provide an estimate of the magnitude of the 
future long term fluctuations which can be expected at the geographical locations cited. The dis- 
tribution of the yearly and diurnal variations is nearly symmetrical and therefore the changes about 
the best yearly average or best daily average of N s will vary approximately as ±An/2, where AN is 
the total expected variation indicated in Table 1. 

It is indicated in the section entitled "Tropospheric Errors” how often refractivity data must 
be obtained to assure attainment of a specified degree of accuracy in elevation angle, range, and 
range-rate measurement. With reference to Table 1, it is seen that a diurnal change in N s of 20 
is seldom exceeded for the sites listed and hence a variation of ±10 or less about the best daily 
average can be anticipated. Similarly, a yearly change of 56 in N s is seldom exceeded and a 
variation of ±28 or less about the best yearly average can be expected for the seven geographical 
locations listed in Table 1. The maximum error remaining after a correction is made for tropo- 
spheric effects is examined by means of ray-tracing techniques for the following cases: 

1. A worldwide "standard” correction of N s = 313. Error in knowledge of N s taken as ±60. 

2. A best yearly average of N s at a particular site. For this case the error in knowledge of 
N s is assumed as ±28 from the actual value at any given time. 

3. A best daily average of N s at a particular site. For this case the error in knowledge of N s 
is assumed as ±10 from the actual value at any given time. 

At this writing, numerous studies are being undertaken in an effort to simplify the task of ob- 
taining meaningful bias correction terms for application to tracking data. For example, the Air 
Force Cambridge Research Laboratory, The National Bureau of Standards, and the Mitre Corpora- 
tion are currently conducting experimental programs in an attempt to attain improved correlation 
between atmospheric structure and atmospheric refraction effects. Cambridge Research Labora- 
tory is trying to arrive at a system whereby the refraction profile of the troposphere can be more 
accurately determined by observing prevailing cloud conditions. This laboratory has put meteoro- 
logical equipment, including a refractometer, aboard a C-130 aircraft which is flown through clouds 
at Cape Kennedy at altitudes up to 6 km. At the same time a U-2 aircraft flies above the C-130 
and photographs the C-130 as well as the prevailing cloud cover. Later the cloud pattern is matched 
with refractometer readings. It is hoped that a sufficient number of comparisons will be obtained 
to establish an overall refractivity profile which includes generally neglected effects of cloud cover. 
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Table 1 

Range of Variation in Surface Refractivity (N s ). 


Weather 

Yearly 
Variation 
in N s , AN 



Average Diurnal Variation in N s , AN 




Yearly 

Average, 

N s 

Station 

Jan. 

Feb. 

March 

April 

May 

June 

July 

Aug. 

Sept. 

Oct. 

Nov. 

Dec. 

Miami, 

Florida 

44 

10 

6 

6 

11 

13 

12 

12 

13 

11 

10 

10 

10 

359 

Flagstaff, 

Arizona 

41 

6 

9 

8 

11 

14 

17 

14 

14 

16 

16 

10 

6 

246 

Brownsville , 
Texas 

52 

10 

12 

14 

17 

16 

19 

23 

22 

15 

14 

11 

11 

360 

Tampa, 

Florida 

55 

14 

13 

14 

16 

22 

17 

14 

12 

16 

17 

12 

8 

356 

Washington, 
D. C. 

56 

6 

8 

9 

10 

11 

14 

18 

16 

16 

12 

9 

6 

334 

Fresno, 

California 

39 

6 

11 

14 

20 

24 1 

25 

26 

27 

24 

17 

9 

6 

318 

Bismark, 
North Dakota 

39 

4 

4 

4 

9 

11 

10 

12 

13 

12 

10 

6 

4 

308 


The National Bureau of Standards is carrying out tests in Hawaii between a 3 km mountain 
peak and a point at sea level 24 km distant. These tests are being made to further improve the 
accuracy to which a troposphere refractivity profile based on variables such as temperature, pres- 
sure, wind, and moisture content can be calculated. The Bureau of Standards hopes to compile 
another series of refractivity tables based on the new findings. 

The Mitre Corporation is currently developing a M line integral refractometer ft to provide a 
measure of refractivity based on only the moisture content (largest unknown in terms of radiowave 
refraction) and oxygen content of the troposphere (Reference 8). 

Finally, the section entitled "Ionospheric Errors" includes some current thoughts regarding 
ionospheric effects above 100 Me. It should be pointed out that azimuth angle measurement errors 
due to the atmosphere's departure from spherical stratification are not considered in this report. 
Measurements at Cape Kennedy demonstrate that for elevation angles above 5° there is generally 
a negligible effect on ray-bending due to tropospheric horizontal inhomogeneity (Reference 9). A 
horizontally inhomogeneous ionosphere will occur whenever a measurement ray passes from a 
daylight to nighttime region; however, the effect on azimuth angle measurement at frequencies 
above 1 kMc is negligible. 


TROPOSPHERIC ERRORS 
General 

Figure 4 indicates the physical relationship between the troposphere and ionosphere above the 
earth's surface. The index of refraction, n, in the troposphere is given by 

n = ^ = y e r K . (1) 
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TARGET ABOVE 
IONOSPHERE 



Figure 4— -Relationship between 
troposphere and ionosphere. 



where 

c = speed of light 

v p = velocity of propagation in medium of 
index n 

e r = relative dielectric constant of me- 
dium and 

fM r = relative permeability of medium. 

For all practical purposes the relative perme- 
ability of the troposphere is that of freespace 
(i.e., n r = 1.0) and therefore the dielectric 
constant at any point within the troposphere is 
a measure of the index of refraction at that 
point. This relationship between e r and n is 
sometimes used in measuring tropospheric 
index of refraction profiles by means of a mi- 
crowave refractometer (Reference 9). 

Another scheme for measuring the lower 
troposphere refraction profile is by means of 
an expendable balloon-borne radiosonde which 
transmits meteorological data back to earth 
from heights of 0 to 15 km. The radiosonde is 
automatically tracked by a high- gain antenna, 
thus assuring sufficient data for a continuous 
profile. Measured profiles for Valkaria, Florida, 
as obtained from this system (the tT Rawin Sys- 
tem") are shown in Figure 5. The two curves 
for 1100 hours EST 7 August 1963 and 1100 
hours EST 8 August 1963 indicate the degree of 
variation one can expect in day-to-day measure- 
ments for the same time of day (i.e., no diurnal 
effect). The slight differences between the pro- 
files are attributed to day-to-day changes in 
cloud cover. 


Figure 5— Measured tropospheric refractivity 
(data from Reference 5). 


Since the index of refraction, n, associated with the atmosphere is always near unity, it is 
general practice to describe this index in terms of the "refractivity," n, where 

N = (n - 1 ) 10 6 - (2) 

This notation is useful when considering either the troposphere or the ionosphere. 
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The dielectric constant (and therefore N) of the troposphere is a function of temperature, total 
atmospheric pressure and water vapor pressure. The empirical relationship used for calculating 
N is given by the National Bureau of Standards (Reference 7) as 


N 


77.6 


P + 


4810 e RH ^ 

-r—)' 


(3) 


where 

T = temperature in degrees Kelvin 
P = total atmospheric pressure (millibars) 

RH = % relative humidity 

e s = saturation vapor pressure (millibars). 

The constants in Equation 3 were determined from a consideration of recent microwave and 
optical determinations of the refractive index of air and the result is considered to be accurate to 
within 0.5% in N for frequencies up to 30 kMc in the ranges of temperature, pressure, and humidity 
normally encountered. Thus time delay and elevation angle errors due to tropospheric refraction 
are frequency independent in the 3 Me to 30 kMc radio frequency range. 

An exponential mathematical model for N within the troposphere has been introduced by the 
National Bureau of Standards. It is of the form 

N = N s exp [-k(h-h,)] . ( 4 ) 


where 

N s = refr activity at earth 1 s surface 
k = decay constant (per km) 

h = height above surface corresponding to N (km) 
h s = height of reference point above mean sea level (km). 

The decay constant is generally calculated from a refractivity measurement at a height of 1 km. 
That is, 


k 


ln N 

S 


N 

s 

tan * 


(5) 


N s +AN = refractivity at height of 1 km above surface (Reference 4). Equation 4 has been 
shown to accurately account for the mean tropospheric refraction of solar radio energy at 
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16.2 kMc and 34.5 kMc, frequencies well 
above the range affected by the ionosphere. 
This measured data was presented previously 
in Figure 1. 

The total measured residual elevation bias 
error at BDA, Bermuda, associated with an 
FPS/16 radar (5400-5900 Me) track of the sec- 
ond orbit of Mercury MA7 is presented in 
Figure 6. While the method for obtaining the 
absolute calibration of this particular radar is 
not presented in Reference 10, it is of interest 
to note that there are no large deviations from 
the observed minus computed (o-c) angle er- 
ror curve. At the time these data were taken, 
only yearly average values of surface refrac- 
tivity were in use at the various tracking sta- 
tions and the magnitude of the residual eleva- 
tion error is not surprising. The measured residual of Figure 6 was obtained using the same 
computed orbit for each data point. 
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Figure 6— Total measured residual elevation angle bias 
for FPS/16 track of MA7. Tracking time from 1 .7° eleva- 
tion to 45.2° elevation is approximately 2 minutes. 


Tropospheric Elevation Angle Error 

The elevation angle error due to tropospheric refraction can be deduced from Snell’s law for 
spherically stratified media. This law is given by 


nr cos e - n 0 r 0 cos e Q 


(6) 


where 


n 



e and e 0 


index of refraction at lower layer 
index of refraction at upper layer 
radius to lower boundary associated with n 
radius associated with boundary of n and n 0 
ray local elevation angles. 


The geometry is indicated by Figure 7. Equation 6 is used as the basis for the derivation of the 
classic expression for the bending of a radio ray (References 4 and 9), namely 


-J> 


cot e 


o* 


( 7 ) 
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It can be shown (Reference 9) by integrating Equation 7 by parts that the bending through the total 
troposphere at elevation angles, e (Figure 8), greater than 5°, the total refraction, r, is given to 
within 10% by 


r - (N s cot 6 ) 10~ 6 radians 


( 8 ) 


where 

r = total refraction 
N s = surface refractivity 
e = ground elevation angle 


(10% accuracy for e = 5°; 
2% accuracy for e > 10°). 


The tropospheric refraction correction currently being employed in the Gemini program is an 
approximate evaluation of the integral of Equation 7. This approximation offers good agreement 
with the National Bureau of Standards ray -tracing tables for elevation angles greater than 2° and 
has the added advantage of being easy to implement into a computer program. This formulation 
was first presented by H. E. Clark of the Data Operations Branch, Goddard Space Flight Center, in 
1963 (Reference 10). With reference to Figure 8, 


[(N s x 10" 6 ) (cot e)] 1 


1. 072014 x 10" 
1.03585796 - -- -z 


1.279119 x 10 


1. x 10 


2° < e < 10 c 


(9) 



Figure 7— Ray bending geometry. 


T = (N s X 10“ 6 ) cot € 10° <€<90°, 
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where 


e = observed elevation angle (radians) 

N s = surface refractivity 
r = refraction (radians). 

The present refraction corrections utilizing Equation 9 are based on the monthly average of 
surface refractivity for a particular site. Since the monthly change in surface refractivity for a 
particular time of day is a rather slow function of time (refer to Figure 3), any deviation from the 
monthly average is primarily due to the cyclic diurnal variations which, as shown by Table 1, will 
seldom exceed a change in N of ±10. 

From the geometry of Figure 8, it is seen that the limits of the elevation angle error, Ae , in 
terms of the refraction, r , are given by 


T 

1 


< A e < t . 



ELEVATION ANGLE, e (degrees) 


Figure 9— Elevation angle error due to tropospheric re- 
fraction. Vehicle height is 200 km. Maximum remaining 
error after correction is based upon exponential tropo- 
sphere ray-tracing. Results are applicable to most 
tracking sites and In particular, those near Miami, 
Florida; Flagstaff, Arizona; Brownsville, Texas; Wash- 
ington, D. C.; Fresno, California; Bismark, North Dakota; 
and Hamilton, Bermuda. 


Furthermore, as the height, h, is increased 
without limit, Ae approaches r. For all prac- 
tical purposes, when considering heights above 
200 km and elevation angles above 5°, r = Ae . 

The NBS ray-tracing tables (Reference 4) 
present in tabular form the tropospheric ele- 
vation angle error, ]Ae |, resulting from ray 
bending in an exponential atmosphere. The re- 
sults are tabulated for each of 8 values of N s 
(i.e., N s = 200.0, 252.9, 289.0, 313.0, 344.5, 
377.2, 404.8 and 450.0) with surface- to- target 
heights from 0 to 450 km. Since 450 km is well 
into the ionosphere (see Figure 4) the higher 
altitude calculations (as pointed out in the sec- 
tion entitled ’Ionospheric Errors”) will re- 
quire modification for an estimate of total er- 
ror. This is especially true if frequencies 
below 2 kMc are employed. 

From the foregoing data the following fig- 
ures have been derived: Figure 9 shows the 
maximum error, \rf € \ , in elevation angle de- 
termination as a function of elevation angle for 
corrections based upon the standard atmosphere, 
best yearly average, and best daily average 
values of N s . 
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The curve marked "no correction" represents a ray tracing using the National Bureau of 
Standards exponential model for a typical value of N s (373 in this case). An examination of the ray- 
tracing data (Reference 4) reveals that the refraction error for a given elevation angle is, for all 
practical purposes, a linear function of surface refractivity for the case of elevation angles greater 
than 2 degrees when target heights are greater than 200 km and surface refractivity is greater than 
200. This means that at any given elevation angle an error in the determination of the surface re- 
fractivity will produce a certain error in elevation angle determination independent of the actual 
value of the surface refractivity. That is, the error after correction is only a function of the dif- 
ference between the true value, N s , and assumed value, N s . 

As indicated in the introduction, the maximum diurnal variation is on the order of ±10 for the 
7 sites listed in Table 1. Also, a value of ±28 was shown appropriate for maximum yearly variations 
in N s . If a standard atmosphere of N s = 313 is used, the uncertainty in N s at any given site can 
often be ±60. Figure 9 is therefore appropriate for most tracking sites. The data in Reference 9 
indicate that the variations in yearly N s at Woomera, Australia, and Madrid, Spain, would be less than 
that experienced at the tracking sites in Florida and hence the results are somewhat pessimistic 
for these locations. Figure 9 is appropriate for estimating the tropospheric introduced elevation 
angle error for satellites or space vehicles at or above an altitude of 200 km. Figure 4, which re- 
lates the physical relationship between the troposphere and ionosphere, shows why total tropospheric 
ray bending is independent of altitude above 200 km. The value of N at 200 km is an inverse function 
of frequency squared and for sufficiently high frequencies can be taken as zero (i.e., freespace). 


Tropospheric Range Error 

The tropospheric range measurement error between a target and the earth’s surface is pri- 
marily due to the electrical path length difference between a line-of-sight path measured in terms 
of the velocity of light in vacuum and the line-of-sight electrical path length due to the actual ve- 
locity of propagation associated with the troposphere. A second-order effect is the physical path 
length difference between the ray path and the line-of-sight path. The National Bureau of Standards 
ray tracings include both factors. A certain degree of insight regarding the major source of tropo- 
spheric range error can be obtained from the following derivation which gives a range error ap- 
proximation (97% of total error for e > 10°) due only to electrical path length difference. 

The path length between ground station and target (points 1 and 2, Figure 8) is given by 


R 0 = J ds = line-of-sight range . (10) 

The "apparent range" due to a velocity of propagation less than that of light in vacuum is given by 


R = 


f 


nds - apparent range 


(ii) 


n 



The difference between Equation 10 and 11 is the one-way range error, AR 7 or 


AR 


f 


(n - 1) ds ; 


( 12 ) 


but 


(n-l) s N(l0" 6 ) 


and Equation 12 becomes 


AR = 



But from the geometry of Figure 7 it is seen that 


ds 


dh 

sin e 


dh esc € . 


(13) 


(14) 


Combining Equations 13 and 14 


AR 


10' 6 


r 


N esc e dh . 


(15) 



Figure 1 0— Tropospheric range error remaining after cor- 
rection versus elevation angle. The error remaining after 
correction represents the difference between ray-traces 
at N s = 377 and N s = 313 (standard atmosphere). 


Equation 15 indicates that the range error, like 
the elevation angle error, is a function of the 
refractivity profile. It can be shown (Refer- 
ence 9) that Equation 15 can, for elevation 
angles e > 10°, be approximated by 



(16) 


where 

h L = height of surface above sea level, 
h 2 = target height, 

which for e > 10° accounts for 97% of the range 
error. Figure 10 indicates the range error re- 
maining after correction with a correction 
based only upon the often accepted standard 
value of N s = 313.0. It is apparent that for 
elevation angles greater than 5°, the range 
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correction is not a critical function of N s since, as shown by Figure 10, the remaining error after 
correction is only 0.5 meters at 5° even though the ray-tracing is for a difference of ±64 between 
the true and assumed values of N. The insensitivity in range error to surface refractivity when 
the entire troposphere is traversed (h > 85 km) is due to the fact that the absolute value of the decay 
constant, k, of Equation 4 is observed in nature to be inversely related to the surface refractivity, 
N s . Thus the integrated time delay experienced in tracing a ray through the entire troposphere is 
to a certain extent invariant with changes in surface refractivity. 

The tropospheric error indicated is appropriate for lunar distances as well as earth orbits. It 
should be pointed out that at lunar distances the limit on range accuracy is the accuracy to which 
the velocity of light is known. An estimate of the range and angle non- systematic or noise errors 
due to tropospheric anomalies in cloud cover is indicated in Table 2 (Reference 12). The noise 
errors indicated in Table 2 are pessimistic in that short term fluctuations are smoothed by the 
overall time response of the tracking system. That is, the system discerns the mean value of N s as 
indicated by the measured data of Figure 6. 


Table 2 

Estimated Range and Angle Noise Due to Troposphere. 




5 Degree 

Elevation 

20 Degree Elevation 

Type of 
Weather 

Refractivity 
Fluctuation, 
AN RMS 

Elevation 
Fluctuation, 
S € RMS 
(mr) 

Range 

Fluctuation, 
S R RMS 
(meters) 

Elevation 
Fluctuation, 
S e RMS 
(mr) 

Range 
Fluctuation, 
S R RMS 
(meters) 

Heavy Cumulus 

30 

0.4 

0.5 

0.2 

0.25 

Scattered Cumulus 

10 

0.2 

0.12 

0.1 

0.06 

Small 

Scattered Cumulus 

3 

0.1 

0.02 

0.05 

0.01 

Clear Moist Air 

1 

0.04 

0.005 

0.02 

0.0025 

Clear Normal Air 

0.3 

0.02 

0.001 

0.01 

0.0005 

Clear Dry Air 

0.1 

0.01 

0.0002 

0.005 

0.0001 


Tropospheric Range-Rate Error 

The errors in range-rate measurement due to signal passage through the troposphere can be 
attributed to the following causes: 

1. The angular difference between the refracted ray at the vehicle and the ground to vehicle 
line-of-sight. 

2. Incorrect assumptions regarding the magnitude of the refractive index at the space 
vehicle. 
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3. Changes in index of refraction profile during the Doppler frequency measurement interval 
introducing an error when integrated Doppler is averaged over a finite time (usually one 
second or greater). 

4. Short term fluctuations in the index of refraction associated with the vehicle-to-ground 
radio path. Such TT noise” fluctuations include those due to clouds, rain, and changes in prop- 
agation path as the vehicle passes through various strata at a constant radio distance from 
the ground station. 


Since this paper is concerned with vehicles at altitudes above 200 km, the vehicle will be well 
above the troposphere (see for example, Figure 4) and consequently, the first two error sources, 
angular error and local vehicle refractive index errors, are attributed to the ionosphere 
and are discussed under "Ionospheric Range-Rate Error." The tropospheric range-rate errors 
generally considered in the literature are those assuming "instantaneous Doppler frequency" meas- 
urements in which case range-rate error for a target traversing the atmosphere is considered as 
pnly the localized error at the specific point in the medium associated with the space vehicle. 

From this viewpoint, the troposphere acts merely as a refraction medium (Reference 13). In 
practice, however, the Doppler shift associated with a satellite pass of space vehicle motion is 



ELEVATION ANGLE, c (degrees) 

Figure 1 1— * Range -rate error due to troposphere (alti- 
tude, 225 km; circular earth overhead orbit vehicle 
speed, 7. 85 km/sec; error decreases for higher altitudes). 


obtained by counting Doppler cycles of phase 
change over a finite time interval (often 1 sec- 
ond) and obtaining the average Doppler rate for 
this interval. With this practical consideration 
in mind, it is noted that the range- rate error 
due to a different index of refraction profile at 
the beginning and end of the sampling period 
will be equal to the difference between the 
range errors associated with the beginning and 
end of the measurement interval divided by the 
time interval of the measurement. 

Figure 11 shows the remaining range-rate 
error, \vk\, as a function of elevation angle 
after a standard atmosphere (N s = 313) cor- 
rection, after a correction based on a best 
yearly average, and after correction based on 
the best daily average of N s . This curve is for 
a 200 km overhead earth parking orbit. The 
tropospheric range-rate error was derived from 
NBS ray-tracing data (Reference 4) in the fol- 
lowing manner. 

The tropospheric range error versus ele- 
vation angle curves for N s = 377, N s = 367, 
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N s = 349, and N s = 313 were graphically differentiated to obtain curves of the rate of change of 
range error with respect to elevation angle, d(AR)/de, versus elevation angle, e. 

Also a curve of earth central angled, (Figure 8) versus elevation angle, e, was obtained from 
the ray tracing data of Reference 4. It was noted that in the range 1°< e < 12° e is linearly approxi- 
mated by 


e - -1.550 + 0.46 radians , (17) 

where 

e = elevation angle (radians) 

0 = central angle (radians). 

And since the vehicle speed is 7.85 km/sec (nominal Apollo earth parking orbit speed), the angular 
velocity d0/dt , is given by 


d0 

dt 


V 

r 


1. 19 x 10~ 3 radians/second 


(18) 


where 

V = vehicle speed = 7.85 km/sec 

r = earth radius plus altitude of vehicle = 6378 + 225 = 6603 km. 

Combining Equations 17 and 18 it is noted that the elevation angular velocity is simply 


de 

dt 


d0 

(-1.55) at- 


1.85 x 10 3 radians/second 


The change in range per unit time can then be obtained as 


1° < e < 12° . 


ARj - ar 2 

AT 


(^37^) Olt) = (^aP)( -1 - 85 x 10' 3 radians/second) 


(19) 


( 20 ) 


where 


de/dt 

A *l 

ar 2 

AT 

d(AR)/de 


elevation angle rate of change (radians/second) 
range delay error at e x (meters) 
range delay error at e 2 (meters) 
time period of integration (seconds) 

rate of change of range error with respect to elevation angle obtained from graphical 
differentiation of ray -tracing range error versus elevation angle curve. 
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Equation 20 is plotted for N s = 377 in Figure 11. For every elevation angle, the change in 
range error in cm over a one second interval represents the range-rate error in cm/sec. It is 
seen that this range-rate error is independent of the averaging time, At, as long as the slope of 
range error versus elevation angle, d(AR)/de, does not change appreciably during the measurement 
interval. This condition is met for a sampling rate AT up to minutes depending on vehicle angular 
velocity relative to the tracking station. 

In the Apollo Mission the lunar parking orbit is to be maintained at a nominal height of 185 km 
above the moon’s surface at a speed of 2.28 km/sec (Reference 14). The range-rate error versus 
elevation angle curve shown in Figure 11 is not appropriate for lunar orbits since it is based on 
earth orbit geometry and velocity. With reference to Figure 12 it is seen that the angular rate of 
change of the vehicle relative to earth central angle 4 > (assuming a stationary earth) is quite slow 
and is given approximately by 

d0 ^ V _ , 

^ - 6 x 10 6 radians/second , (21) 


where 

V = lunar orbit linear speed = 2.28 km/sec 
R 0 = earth-moon mean distance = 3.84 X 10 5 



Figure 12— Lunar orbit geometry (Apollo lunar parking 
data from Reference 14). 


km. 

Therefore at lunar distance and beyond, 
the maximum range-rate error is due solely to 
the earth’s rotation, which is on the order of 
7 X 10“ 5 radian per second. The uncorrected 
range-rate error at lunar distances for N s = 377 
is shown in Figure 13. 

A verification of the foregoing analysis is 
presented by Figure 14, which shows the esti- 
mated residual two-way Doppler at 960 Me from 
Ranger VI (Reference 15) at a distance where 
only earth rotation is important. The range - 
rate error corresponding to the two-way Dop- 
pler is given by 

dR _ d <j> (k\ _ e (k\ _ 

dt - dt^J " f d \ 2 ) " 15.6 f d cm/sec , (22) 

where 

f d = two-way Doppler rate (cps). 

The calculated range-rate error of Fig- 
ure 13 due to the earth’s rotation is given by 
Equation 20 for the case de/dt = - 7 X 10“ 5 
radian/second. 
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Figure 1 3— Range-rate error due to 
troposphere and earth rotation. 


Figure 14— JPL Doppler correction compared to cor- 
rection based upon ray-tracing (elevation rate ~ 0.245 
deg/min = 7 x 10~ 5 radians/sec). 


It is seen that Figure 13, when replotted in Figure 14, fits the Jet Propulsion Laboratory data 
quite well. The Jet Propulsion Laboratory uses an empirical correction for the tropospheric re- 
fraction effect upon range- rate data. The JPL data in Figure 14 were obtained by generating best 
fit orbit data and then comparing these data to those obtained without the refraction correction 
terms. It appears that the ray-tracing range- rate refraction correction derived in this paper is 
appropriate for sampling times at least up to one minute, the sampling rate associated with Ranger 
VI in Figure 14. It is seen that rather straightforward correction terms can be obtained from the 
derivatives of the ray-traced range error curves for either earth orbit or lunar distance space 
travel. At the particular frequency of 960 Me used on Ranger VI, the total range- rate error due to 
the atmosphere can be accounted for by the 

troposphere. Table 3 


An estimate of tropospheric short term 
Doppler error (noise) based on range fluctua- 
tions measurements by NBS is presented in 
Table 3 (Reference 16). 

The noise estimate is seen to be consid- 
erably less than the biasing error due to trop- 
ospheric effects and is generally independent 
of elevation angle. The reason for this is il- 
lustrated by the equation describing the Dop- 
pler signal as received at the ground sta- 
tion. Under 1 ’Tropospheric Range Error/ 1 it 


RMS Range- Rate Noise Due to Tropospheric 
Fluctuations (1-Way Transmission). 


Doppler 

Integration 

Time 

(seconds) 

Cloud Coverage 

Small 

Average 

Large 

R Noise (cm/sec) 

1.0 

0.016 

0.022 

0.075 

5.0 

0.015 

0.020 

0.055 

20.0 

0.010 

0.015 

0.030 

60.0 

0.003 

0.009 

0.020 
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was shown that the apparent range, R , is given by 


R = 


f 


n ds 


p<t) 

I n(t , s) ds 

J Sl (t) 


(23) 


where 


n(t , s) = tropospheric index of refraction along path of integration, a function of time and path 
s 2 = ray terminating point at vehicle 

Sj = ray terminating point at earth = 0 or reference of coordinate system used in per- 
forming the integration and the geometry is as indicated by Figure 8. 


But the range rate dk/dt = R is given by 


d p (t) f* /3n\ ds l 

Tt n(t, s) ds - \3t / ds ~ n i ~dF + n 

J s j(t) = 0 JO 


2 '3n\ . ds i ds 2 

2 dt 


or 


= p® 

Jo 


ds + n„ 


ds 2 
dt 


(24) 


Thel first term in Equation 24 relates the measured "apparent change in range rate" due to the 
fluctuations in index of refraction, n , along the ray path between points 1 and 2. The short term 
fluctuations are generally due to cloud cover and hence not a critical function of elevation angle e . 
The first term includes the change in n during Doppler averaging. The second term is a function 
of the range rate dsy'dt, which is the velocity of the vehicle projected upon the tangent of the ray 
path at the vehicle, and the index of refraction at the vehicle. The second term is entirely a func- 
tion of the medium characteristics at the vehicle and for a vehicle height of 200 km is concerned 
with the ionosphere. 


IONOSPHERIC ERRORS 


General 

Correction procedures based on use of ionospheric profiles have not proved very accurate to 
date (1965) and may not lead to any significant reduction of error. The estimates of error reduc- 
tion run from negligible to a factor of four (i.e., leaving 25% of the initial error). The improvement 
realized is quite dependent upon the stability of that portion of the ionosphere through which the 
measurement ray passes. The ionospheric layers are always in a state of motion near dusk or 
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dawn (Reference 12). It is therefore desirable to obtain an estimate of the maximum daytime or 
nighttime ionospheric error which could be experienced at any given frequency and consider this 
value as the limiting uncorrectable error due to ionospheric effects. 

The index of refraction at a given point within ionized media such as the ionosphere is defined 
as the ratio of the velocity of light in vacuum, c , to the phase velocity, V p , or 


n 


c 


V 


p 



N e 2 


1/2 


(Reference 17) 


(25) 


where 

V p = phase velocity 
c = speed of light in vacuum 
N e = electron density = electrons/meter 3 
e = electron charge = 1.602 x 10" 19 coulombs 
m = electron mass = 9.11 X 10“ 31 kilograms 

e 0 = free space dielectric constant = 8.855 x 10~ 12 farad/meter. 

This definition of index of refraction is consistent with that presented in Equation 1 for non- 
conducting media where "velocity of propagation" is equal to "phase velocity." In order to properly 
interpret Equation 25 the concept of phase and group velocity must be briefly reviewed. 

Phase velocity is that velocity one must travel to keep the instantaneous phase constant; that is, 
if e j(&,t ~^ z) represents the phase factor associated with a sinusoidal continuous signal where 

oj = angular frequency = 2tt f 
f = frequency of oscillation 
J3 = phase delay per unit distance = <^/V p 
V p = phase velocity 

z = distance between source and observer, 

then it is noted that the instantaneous phase t - f3z) does remain constant if one travels with 
velocity V p such that z = V p t +K, where K = a constant. 

In nondispersive media (i.e., where the phase velocity is not a function of frequency) such as 
the troposphere the "phase velocity" and "velocity of propagation" (or more precisely signaling 
velocity) are identical. This is not so in the ionosphere as might be deduced from Equation 25 
which indicates | V | > c . 

Signal velocity is associated with the velocity of propagation of the envelope of a modulated 
sinusoid.. The modulation in the case of a ranging system is often simply a pulse or series of 
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pulses. Since a modulated sinusoid will consist of a spectrum of frequencies, phase distortion is 
inherent in dispersive media. However, if the transmitted spectrum is narrow compared with the 
carrier frequency, an average tT signal or group velocity” can be calculated which can be used to 
calculate time delays through the ionosphere. This velocity is given by (Reference 18) 


V 


g 



(26) 


Differentiating Equation 25 and substitution into Equation 26 yield, after some manipulation, 


V 


K 


V n 2 - nc 

p 


(27) 


for 


n < 1 . 

That is, since n < 1, the group velocity in the ionosphere is less than the velocity of light in vacuum 
and hence a time delay in signaling is realized. 


Ionospheric Elevation Angle Error 


In principle, the ray-bending or refraction mechanism within the ionosphere is analogous to 
that of the troposphere. The measurement ray traverses a region of varying refractivity which is 

a function of electron density and excitation 
frequency as indicated by Equation 25. For 
this reason, ray-tracing techniques can be 
employed to obtain estimates of ionospheric 
refraction error. Figure 15 indicates the 
nature of the combined troposphere and iono- 
sphere refraction geometry. The details of 
the overall ray path are related to the overall 
refractivity profile. At present only the tropo- 
sphere refractivity profile is known with suf- 
ficient accuracy to permit real time bias ad- 
justments of radar data. Models based on 
average daytime and nighttime ionosphere 
electron density profiles, however, are useful 

c . ic t l • j in estimating average biasing effects which 

Figure 1 5— Tropospheric and o o & 

ionospheric ray refraction. can be attributed to the ionosphere. Separate 
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ionospheric models are generally used for day and night since daytime ionization is considerably 
greater than nighttime ionization. 

The value of refractivity, N, for the ionosphere is given by (Reference 16) 

- 40.3N e 

T 2 




where 

N e = electron density (electrons/meter 3 ) 
f - frequency (cycles per second) 

and 


N = (n ~ 1) 10 6 . 

An example of a typical ionospheric refractivity profile at 2 kMc based upon data measured at 
Washington, D. C., is shown in Figure 16. These curves were originally presented in Reference 19 
for 400 Me and have been adjusted to 2 kMc by means of Equation 28. The maximum value of 
N = - 12 corresponds to the maximum value of electron density N e = 1.2 x 10 12 electrons per cubic 
meter at an altitude of 300 km. A value of 10 12 electrons/meter 3 for N e is seldom exceeded in the 
earth T s ionosphere; hence a value of N = - 12 is the maximum one might expect during any iono- 
spheric ray tracing at 2 kMc. This is in contrast to the values of tropospheric N at the earth T s 
surface which range from 300 to 400. The course and fine structure of the ionosphere is a subject 
of continuing measurement. Techniques for such investigation by means of satellite signal observa- 
tion are presented in References 20 and 21. 

As in the case of tropospheric refraction 
the maximum error in elevation angle occurs 
at low elevations. Also, maximum refraction 
is associated with maximum ionization and 
hence daytime ionospheric states. Figure 17, 
which is for 5° elevation and a daytime average 
ionosphere, thus relates the maximum mean 
elevation bias error as a function of frequency 
and altitude above the earth T s surface. This 
analysis is based on a rectangular model of the 
ionosphere and was first presented by Pfister 
and Keneshea (Reference 22). Figure 18 indi- 
cates the rectangular electron density versus al- 
titude approximation used to obtain ionospheric 
error estimates. While such a model is useful in 
estimating maximum error, it is apparent that it 
cannot be used for purposes of radar calibration. 


Q 


1200 
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800 
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Figure 1 6— Ionospheric refractivity profiles at 2 kMc. 
Curves are reproduced from Reference 19 with a scale 
change according toEquation28 to indicate the range of 
refractivity at 2 kMc. Profile A exceeded 97.5% of the 
time, B, 50% of the time and C, 2.5% of the time. 
Maximum electron density = 10 12 electrons/m 3 . 
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Figure 17 — Ionospheric maximum elevation angle error 
for 5° elevation (from Reference 22). Error decreases as 
elevation angle increases above 5°. 



Figure 18 — Daytime ionospheric model. 


A number of measurements carried out 
on the California coast indicated that the 
ionospheric errors could not be predicted on 
to obtain current profile data) any better than 


the basis of ionograms (ionospheric soundings 
by using average monthly forecasts of electron density profiles and that the residual errors 
after attempted correction were almost as great as .the original error values. The inability 
to predict the error encountered at a given time indicates that accurate systems must oper- 
ate at frequencies high enough to reduce the initial values of ionospheric error to a tolerable 
value (Reference 1). 


Ionospheric Range Error 

The range bias due to increased time delay relative to freespace delay of the signal through 
the ionosphere can be calculated only if the ionospheric index of refraction profile is accurately 
known. If only an average profile is known, the residual error can be estimated. The delay in an 
ionized media must be calculated using the group or signaling velocity indicated on page 20. The 
ionospheric range error due to time delay can then be calculated in a way analogous to the tropo- 
spheric development on pages 11 and 12; 


R 


o 



1 ine-of -sight range. 


(29) 


The "apparent range" due to a phase velocity greater than that of light in vacuum is given by 


R 



measured range, 


(30) 


where 


n = ionosphere index of refraction. 
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Equation 30 describes the apparent increase in range due to a group velocity less than the speed of 
light. Combining Equations 29 and 30, range error AR is given by 


AR = R - R 0 



but n = 1 + N x 10" 6 and Equation 31 can be written as 


(31) 


“ = " 10 '* { * " l °‘" j, Nds ' (32) 

where, for the ionosphere, 



N e = electron density (electrons/ meter 3 ) , 
f = frequency (cps). 


Figure 19, which is for 5° elevation and a 
daytime average ionosphere (a worst case sit- 
uation), shows the maximum mean range bias 
error as a function of frequency and altitude. 

Figure 20 shows the maximum ionospheric 
range error which can be expected at the Apollo 



AR 

Figure 19— Ionospheric maximum range error for5°eleva- 
tion (from Reference 22). Error decreases as elevation 
angle increases above 5°. 



Figure 20— Ionospheric maximum range 
error at 2 kMc (from Reference 22). 
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RANGE FLUCTUATION (meters) 


Figure 21 —Ionospheric range fluctuation versus 
frequency for a 10° elevation angle. 

Ionospheric Range-Rate Error 


Unified S-band frequencies (nominal 2 kMc) 
versus elevation angle. This figure is de- 
rived from the rectangular ionospheric re- 
fractivity profile analysis presented in Ref- 
erences 16 and 22. 

Figure 21 shows an estimate of iono- 
spheric range fluctuation or noise as a 
function of frequency for a vehicle height 
of 240 km and 10° elevation (Reference 16). 


It was shown on page 18 that the range-rate error can be separated into two parts, that is, the 
error due to time variations of the index of refraction, n, along the propagation path during the 
Doppler measurement interval and the localized error at the vehicle. In the case of a target 
traversing the ionosphere only the localized error is significant in the calculation of ionospheric 
error contribution (Reference 13). The value of the ionospheric index of refraction at the vehicle, 
n t , can only be estimated and presents one limitation on overall Doppler error determination. The 
major ionospheric error source is due to the slight difference in direction of the refracted ray at 
the target and the ground to target line- of- sight. 


The geometry of Figure 8 indicates a target velocity, V, the radial component along the line- 
of-sight, R, and the apparent or measured radial component, V m . The magnitudes of the Doppler of 
R and V m are given by 

R = V cos (0 + 8 ) (33) 

v m = V cos 9 . ( 34 ) 

The magnitude of the error in measuring R is then given by 

AV r = R - V, = V [cos (0+8)- cose] . ( 35 ) 

But cos (<9 + S) = cos 6 cos 8 - sin 6 sin 8 and the angle S is in all practical situations quite small 
such that Equation 35 can be approximated by 

AV r - - VS sin (9 . ( 35 ) 

That is, the magnitude of the range- rate error is proportional to the target velocity, V , the angle be- 
tween the refracted ray at the target and line-of- sight, 8 , and the sine of the local elevation angle, 6 . 
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Figure 22 shows an estimate of ionospheric range-rate bias as obtained by considering 
an average ionospheric profile and Equation 36. The radial velocity error is for Apollo earth 
orbit speeds (7.85 km/sec). The calculation is based on a height of 300 km which is the re- 
gion of maximum ionization (hence maximum bending), as indicated by Figure 16. Therefore 
Figure 22 can be used to estimate the maximum ionospheric residual range-rate bias for 
earth orbiting vehicles traveling at a speed of approximately 8 km/sec. At lunar distances 
the lack of lunar atmosphere results in 8 - 0 in Equation 12 and hence no significant error is 
contributed by the ionosphere. 

Figure 23 depicts the maximum range-rate error attributable to the ionosphere as a function 
of elevation angle for the Nominal Unified S-band System (USBS) frequency of 2 kMc. 


The range- rate bias due to a change in ionospheric range error over the range- rate measure- 
ment interval can be shown to be small relative to the same effect attributed to the troposphere. 
The reason for this is that the range delay experienced by a signal traversing the ionosphere 
is not a critical function of elevation angle. For example, at 1 kMc, at a height of 200 km 
or greater, and for a daytime ionosphere, the maximum rate of change of range error versus 
elevation angle is approximately 6 meters per radian, nearly independent of elevation angle (Ref- 
erence 22). For a 225 km earth orbit and a velocity of 7.85 km/sec, this corresponds to a bias of 
lcm per second at elevation angles between 4° and 12°. It is noted that the tropospheric bias at 


5° elevation due to the same phenomena is 50 cm 
data for vehicles within the ionosphere ever need 
apparent this component of error is not negligi- 
ble. At a frequency of 1 kMc or higher the tropo- 
spheric range- rate error will predominate by a 
factor of 5 or more for elevation angles up to 
15°. In any case, at 1 kMc or above, this par- 
ticular ionospheric range- rate biasing effect 
is less than that which can be attributed to the 
troposphere. 



MAXIMUM RANGE -RATE ERROR (cm/sec) 


Figure 22 — Ionospheric maximum range-rate error (from 
Reference 13) daytime ionosphere; error due to refraction 
at vehicle. 


per second (Figure 11). However, if range- rate 
to be corrected to within 1 cm per second, it is 



Figure 23— Ionospheric maximum range-rate error at 
2 kMc; overhead earth parking orbit, h = 300 km, error 
due to refraction at vehicle (from Reference 13). 
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CONCLUSIONS 


The following is a two-part summary of this report. The first part consists of general state- 
ments concerning the influence of the atmosphere on microwave propagation between the earth's 
surface and a space vehicle located at or above an altitude of 200 km. This portion of the dis- 
cussion examines separately the effect of the atmosphere upon measurements of elevation angle, 
range and range rate. The second part of this discussion summarizes the influence of the atmos- 
phere upon radiowave propagation again for heights greater than or equal to 200 km but at the 
particular nominal radio frequency of 2 kMc. This frequency is of interest in that it approximates 
the operating frequency employed in the Unified S-Band System of the Apollo Program. 

Atmospheric Induced Error for Vehicle Heights Above 200 km 

Elevation Angle 

Elevation angle bias due to the atmosphere decreases with increasing elevation angle. Tropo- 
spheric refraction is generally frequency independent whereas ionospheric refraction is inversely 
proportional to frequency squared. Maximum uncertainty in angle of arrival determination occurs 
for vehicles at low elevations at points outside of the earth’s atmosphere. With no attempt at 
correction, the maximum elevation angle bias will be on the order of 4 mr at 5° elevation inde- 
pendent of frequency for frequencies above 1 kMc (Figures 9 and 17). With a tropospheric re- 
fraction correction based on the best daily average of surface refractivity, the tropospheric bias 
at 5° elevation will be on the order of 0.1 mr (Figure 9). A 0.1 mr bias at 5° elevation attributable 
to the ionosphere can occur during the daytime at a frequency at or below 900 Me (Figure 17). The 
nighttime maximum bias effect due to the ionosphere is about one third that experienced at daytime. 

Elevation angle short-term fluctuations (period of fluctuation up to 10 minutes) also decrease 
with increasing elevation angle and for the troposphere are a function primarily of cloud cover. 

At frequencies above 1 kMc, the short-term fluctuations are associated primarily with the tropo- 
sphere and estimates range from 0.4 mr RMS for heavy cumulus cover to 0.01 mr RMS for clear 
dry air. At 20° elevations, the corresponding estimates are 0.2 mr RMS and 0.005 mr RMS 
(Table 2). 


Range 

Range biasing effects due to the atmosphere decrease with increasing elevation angle. The 
tropospheric range error is frequency independent while ionospheric range error is inversely 
proportional to frequency squared. With no attempt at correction the tropospheric range bias is 
approximately 24 meters at 5° elevation. With a correction based only upon the worldwide stand- 
ard atmosphere (N s = 313.0), this range error can be reduced to less than 0.5 meter at 5° ele- 
vation (Figure 10) for any of the typical tracking areas cited in Table 1. For a daytime ionosphere 
a maximum range error of 0.5 meter RMS can occur at 5° elevation at frequencies below 3.5 kMc 
for a 240 km earth orbit and below 10 kMc for vehicle heights at or above 500 km. For a nighttime 
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ionosphere a maximum range error of 0.5 meters RMS can occur at 5° elevation at fre- 
quencies below 2.0 kMc for a 240 km earth orbit and below 5.0 kMc for vehicle heights at 
or above 500 km (Figure 19). At lunar distances the range bias is determined by the un- 
certainty in the velocity of light of one part in 10 6 which results in a bias on the order of 
380 meters. 

Short-term range fluctuations (period of fluctuation up to 10 minutes) decrease with 
increasing elevation angle and, for the troposphere, can be attributed primarily to the effects 
of cloud cover. Estimates of tropospheric range fluctuations at 5° elevation range from 0.5 
meters RMS for heavy cumulus cover to 0.02 cm RMS for clear dry air. At 20° elevation 
the corresponding estimates are 0.25 meter RMS and 0.01 cm RMS (Table 2). Above 1 kMc 
the range fluctuation error due to the ionosphere at an elevation of 10° would not be expected 
to exceed 0.2 meter RMS even under worst case daytime conditions (Figure 21). 


Range Rate 

Range-rate bias due to the atmosphere decreases with increasing elevation angle. The 
troposphere introduces biasing for both earth orbit and lunar distance (or beyond) space travel. 
Above 1 kMc the ionosphere need only be considered for the case of earth orbit. In earth 
orbit there are two major sources of range-rate bias. The first has to do with the path 
change through the troposphere during the range-rate measurement interval and leads to a 
range-rate bias on the order of 50 cm per second at 5° elevation. If a correction is made 
based on the best yearly average of N s this can be reduced to 5 cm per second. A correc- 
tion based on the best daily average reduces the error to essentially zero (Figure 11). The 
second range-rate bias source for an earth orbit through the ionosphere is due to the differ- 
ence in angle between the refracted ray and the line-of-sight as measured at the vehicle. At 
1 kMc the maximum error will be on the order of 100 cm per second at 0° elevation and 
60 cm per second at 15° elevation. These values are reduced by a factor of approximately 
three for nighttime operation. 

For space flight beyond the ionosphere, the error in local ray angle is not significant and 
the primary error source is the troposphere. At lunar distances the earth’s rotation of 
7 x 10" 5 radian per second introduces a tropospheric bias on the order of 2 cm per second 
at 5° elevation (Figures 13 and 14) independent of operating frequency. 

Influence of the Atmosphere at the 
Nominal Apollo Unified S-Band Frequency of 2 kMc 

Tables 4 through 7 indicate the remaining error which can be expected after correction 
at a frequency of 2 kMc for day or night operation for both earth orbit and lunar distances. 
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Table 4 

Expected Error in Measurement Due to Propagation Through the Earth's Atmosphere at 2 kMc.* 



RMS Elevation Angle (mr) 

RMS Range (meters) 

RMS Range- Rate (cm/ sec) 

Type of 
Correction 


Noise 


Noise 


Noise 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Average 

Cloud 

Cover 

Heavy 

Cloud 

Cover 

None 

4 

0.02 

0.4 

25 

0.01 

0.5 

75 

0.022 

0.075 

Yearly Average 

0.3 

0.02 

0.4 

2.3 

0.01 

0.5 

30 

0.022 

0.075 

Daily Average 

0.1 

0.02 

0.4 

1.8 

0.01 

0.5 

25 

0.022 

0.075 


♦Daytime ionosphere 
Earth orbit h = 200 km 
5 ° elevation 

ft noise for 1 sec integration 


Table 5 


Expected Error in Measurement Due to Propagation Through the Earth 1 s Atmosphere at 2 kMc.* 



RMS Elevation Angle (mr) 

RMS Range (meters) 

RMS Range- Rate (cm/ sec) 

Type of 
Correction 


Noise 


Noise 


Noise 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Average 

Cloud 

Cover 

Heavy 

Cloud 

Cover 

None 

4 

0.02 

0.4 

24 

0.004 

0.5 

58 

0.022 

0.075 

Yearly Average 

0.3 

0.02 

0.4 

1.0 

0.004 

0.5 

13 

0.022 

0.075 

Daily Average 

♦Nighttime ionosphere 

0.1 

0.02 

0.4 

0.5 

0.004 

0.5 

8 

0.022 

0.075 


Earth orbit h - 200 km 
5 ° elevation 

ft noise for 1 sec integration 

Table 6 


Expected Error in Measurement Due to Propagation Through the Earth T s Atmosphere at 2 kMc.* 



RMS Elevation Angle (mr) 

RMS Range (meters) 

RMS Range- Rate (cm/ sec) 

Type of 
Correction 


Noise 

■ 

Noise 


Noise 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Average 

Cloud 

Cover 

Heavy 

Cloud 

Cover 

None 

4 

0.02 

0.4 

34 

0.01 

0.5 

2 

0.022 

0.075 

Yearly Average 

0.3 

0.02 

0.4 

10.5 

0.01 

0.5 

0.2 

0.022 

0.075 

Daily Average 

♦Daytime ionosphere 

0.1 

0.02 

0.4 

10 

0.01 

0.5 

Bias 

<0.2 

0.022 

0.075 


Lunar orbit 
5 ° elevation 

ft noise for 1 sec integration 


28 



Table 7 


Expected Error in Measurement Due to Propagation Through the Earth's Atmosphere at 2 kMc.* 



RMS Elevation Angle (mr) 

RMS Range (meters) 

RMS Range- Rate (cm/ sec) 

Type of 
Correction 


Noise 


Noise 


Noise 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Clear 

Normal 

Air 

Heavy 

Cumulus 

Bias 

Average 

Cloud 

Cover 

Heavy 

Cloud 

Cover 

None 

4 

0.02 

0.4 

27 

0.004 

0.5 

2 

0.022 

0.075 

Yearly Average 

0.3 

0.02 

0.4 

3.5 

0.004 

0.5 

0.2 

0.022 

0.075 

Daily Average 

0.1 

0.02 

0.4 

3 

0.004 

0.5 

Bias 

<0.2 

0.022 

0.075 


♦Nighttime ionosphere 
Lunar orbit 
5° elevation 

ft noise for 1 sec integration 
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Appendix A 

Glossary 


For purposes of this paper the following definitions are employed: 
REFRACTIVITY 

The refractivity of the atmosphere is represented by N, where 

N = (n - 1) 10 6 , 
n = index of refraction. 


TROPOSPHERE 

That region of the earth’s atmosphere immediately adjacent to the earth and extending upwards 
about 30 km. 

FREESPACE REGION 

The index of refraction of the earth’s atmosphere between 30 km and 85 km is for all practical 
purposes unity and hence this region can be considered as "free space.” 

IONOSPHERE 

That region of the earth’s atmosphere in which the constituent gases are ionized by radiation 
from outer space. This region extends from about 85 km to 1000 km above the earth. 

ATMOSPHERE 

The gaseous matter surrounding the earth, including both the troposphere and ionosphere. 
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